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Abstract 

In this paper we present the exact solution for the average minimum energy of the random 
bipartite matching model with an arbitrary finite number of elements where random paired 
interactions are described by independent exponential distribution. This solution confirms the 
Parisi conjecture proposed for this model earlier, as well as the result of the replica solution 
of this model in the thermodynamic limit. 
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1 The Model 



The model under consideration can be formulated as follows. We have a society consisting of N 
"men" (labeled by i = 1,2, ...,N) and N "women" (labeled by j = 1,2,...,N) described by a 
given set of N 2 random non- negative interactions {Jij} between every man and every woman. The 
statistics of Jy's is defined by a probability distribution function P[Jij]. 

Then we consider all possible "marriages" with the strict monogamy: every man can be con- 
nected with one and only one women, and vice versa. Thus, a particular marriage configuration in 
this society can be described by the N x N permutation matrix Sij with the elements taking values 
or 1 ("0" for all non-coupled pairs of men and women, and "1" for married couples) constrained 
by two conditions: 



N N 

E^- = E^ = i (i.i) 

i=l j=l 

which allow one and only one "1" in each row and in each column of the matrix S. The total 
number of all possible marriage configurations in this society is thus equal to N\. 

Now for every marriage configuration S we introduce the total energy, or total weight (the 
Hamiltonian): 



N 

H[S;J]= J2 SijJii (1-2) 

For a given matrix S this energy is equal to the sum of N particular J^-'s (one from each line and 
each column) corresponding to the particular married couples. In this paper we consider the simplest 
possible model in which the interactions { J^} are assumed to be independent and described by the 
bounded exponential distribution: 



P[Jij] = II ex P (-%) ; (° < Jij < +°°) ( L3 ) 

The problem studied below is formulated as follows: one has to find the value En of the average 
(over the distribution P[Jij}) minimum (over all configurations of the permutation matrix S^) 
energy QO]): 
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P[^]min( Y. S tJ J l3 \ (1.4) 



Equivalently, in the language of statistical mechanics E^ can be obtained as the zero-temperature 
limit of the average free energy: 
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Y[ / dJijexp(—Jij) 
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-f3 Y^ SijJr 
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(1.5) 

Thus, we face the typical problem of statistical mechanics with quenched disorder: first, for given 
values of random parameters {Jij} one has to compute the partition function and the free energy, 
and only after that one has to carry out the averaging over J^'s. 

In the thermodynamic limit (N — > oo) this problem has been solved some years ago in the 
framework of the replica symmetric ansatz yielding the result: 

fljv-oo = C(2) = %- (1.6) 
o 

In this paper we present the exact solution of this problem for an arbitrary (finite) value of N. 
The case N = 1 is trivial: 




E N= i = 1 (1.7) 

The case N = 2 is only slightly more complicated, and it can also be easily calculated explicitly. 
Here the 2x2 permutation matrix S can have only two configurations: 



and 
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Thus, according to the definitions (|1.4j) or (|1.5|) we have: 



/■oo 

E N=2 = 2 dJ n dJ 12 dJ2idJ 2 2(Jn + J22) exp{— J n - J 12 - J 2 i - ^22)^(^12 + J21 ~ Jn ~ J22) (1-10) 
Jo 

Here the ^-function ensures that the state ( |1.8|) has lower energy than the one ( |1.9| ) (due to obvious 
symmetry of the system the contribution from the opposite situation turns out to be the same, and 
this provides the factor 2 in the above equation). Simple integration yields: 

E N=2 = 1 + i (1.11) 
Noting that the result (|1.6|) for N = 00 can be represented also in the form: 

00 j 

Ejt^o = C(2) = Eu ( L12 ) 
k=i K 

and taking into account the results (|1.7|) and (|1.11|) , G.Parisi has recently proposed very elegant 
conjecture that the solution of the problem for arbitrary value of N must be the following 0: 
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;i.i3) 



For the direct calculation of En (in the style of eq.( |1.10|) ) with an arbitrary N one should 
perform the integration over the parameters {Jij} in the constrained positive subspace J^ > 
of the iV 2 -dimensional space. Since the total number of states of the N x N permutation matrix 
is equal to iV! this integration is also constrained by (N\ — 1) hyperplanes which guarantee that 
chosen one particular state has the minimum energy. One can easily verify that even in the case 
N = 3 such calculation turns out to be extremely difficult problem. Nevertheless, simple numerical 
tests for iV = 3,4, 5 proved to be compatible with the above conjecture with the precision ~ 10 -5 
0. Moreover, recent analytical studies have provided the exact solution of this problem up to 
N = 4 and the result of this solution confirms the conjecture (|1.13|) 0. Here we use the original 



idea (proposed by S.Bravyi) of this unpublished work to prove that the conjecture ( |1.13|) is indeed 
correct for arbitrary N. 



2 The Proof 



To ease further presentation of the proof let us introduce the following notation. The operation of 
finding the average of minimum energy of the N x N problem (defined in eqs. ( |1. 4| ) 
be denoted by the symbol 



or 



L5|)) will 
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E 



N 



(2.1) 



It is assumed that "empty" boxes in the above matrix actually contain random elements {Jij} 

Let us consider the first line of the random matrix J^, and among N its elements J\j let us 
find the minimum one: J^ 1 ' = minj(Jij). Due to obvious symmetry of the problem with respect to 
permutations of the columns of the matrix we can always place this minimum element in the 
position (1, 1). Now let us redefine the elements of the first line as follows: 



Jy = J (1) + J Xj , (j ^ 1) (2.2) 

and leave all the other elements unchanged. According to (|1.3j ), the elements Jij are described by 
the same exponential distribution: P[J\j\ = exp (—Jij) , (Jy > 0), while for J^ the distribution 
is: 



P[jW] = N exp (-iVJ (1) ) (2.3) 

Due to the constrains (|1 . 1| ) the above redefinition produces only simple shift of the Hamiltonian 
: 
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N 



H = J« + £ S l3 J K 



(2-4) 



where, the random matrix J^- contains "0" in the position (1,1), while the rest of its elements are 
described by the same distribution (|1.3|). Now using the definition of E N , eq.(|1.5|), we can easily 
integrate out to get: 



where 



En = N + E ^ 



(2.5) 
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(2.6) 



/ 



To calculate E$ let us consider the second line of the above random matrix, and among iV its 
elements J<ij let us find the minimum one: = minj(J2j). Now, due to "0" in the position (1,1) 
the first column of this matrix is no more equivalent to the rest of (N — 1) columns (which remain 
to be equivalent among themselves). Therefore, with the probability 1/N the minimum element can 
be in the position (2,1), and with the probability (iV — 1)/N it can be in the rest of the positions of 
the second line, and in this last case we can place it in the position (2,2). Then we shit the values of 
the elements of the second line: J23 — + J 23 (which leave the distribution of { J23} unchanged). 
The integration over 

J(2) 

gives one more factor 1/N, and for En we get: 
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N 



. (N - 1) (2) 1 ~ (2) 
E N -- + ~^^E N + -E N 



[2.7) 



where 
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Eq.( [2.7|) can be represented in the form: 
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(2-10) 
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To calculate the value E$ defined by the matrix: 
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(2.11) 






































































(2.12) 



let us consider its third line, and among N elements J$j let us find the minimum one: J( 3 ) EE 
miiij (Jsj). Due to two "0" in the positions (1,1) and (2,2) the first and the second columns of 
this matrix are equivalent between themselves, but they are not equivalent to the rest of (JV — 2) 
columns (which remain to be equivalent among themselves). Therefore, with the probability 2/JV 
the minimum element can be placed in the position (3,2), and with the probability (JV — 2)/JV it 
can be in the rest of the positions of the third line, and here we can place it in the position (3,3). 
Then we shit the values of the elements of the third line: J^j = + J$j (which again leave the 
distribution of {^3^} unchanged), and integrate over J® which gives one more factor 1/JV. In this 
way we get: 



En ~ n + — 



-E 



(3) 



2 ~ 



N + jjE® + 5E 
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(2.13) 



where 
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and 



(2.14) 
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Eq.( |2.13| ) can be represented in the form: 



(2.15) 



P _ 3 , Ei(3) , 2 f P (3) , 1 rp(2) 

N N N N N N 



where 
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Proceeding in this way up to the last line we eventually get: 



(2.16) 



(2.17) 
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k=2 



(note that E$ } = since it is given by the matrix with all zeros on the diagonal) where 



(2.1* 
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(2-19) 

Here the double lines mark the positions of the k-th column and the k-th line. 

It can be proved (see Appendix A) that the above value SEffi is given by the rectangular N x k 
random matrix problem: 



6EP = E 



defined by the Hamiltonian: 







































































































(2.20) 



N k 

»=i j=i 



'j 



(2.21) 



where the random matrix is shown in eq.( |2.20D (with the same independent exponential distri- 
butions of non-zero elements). Here the "truncated" N x k part of the original permutation matrix 
S again can have only one " 1" in each line, and besides it has k columns each containing only one 
"1" and (N — k) columns each containing only "0". 

It turns out that the above "rectangular" problem, eq. (|2.20|) can be solved explicitly (the proof 
see in Appendix B): 



SEP 



k-l 



I 



k(k-l)^N-l 
Substituting this result into eq.( [2.18|) we find: 



(2.22) 



z, , 1 £ 1 ^ 1 
EN = l+ N^ 2 k^—l 



(2.23) 



After simple algebra one eventually finds: 
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E N - E N ^ = ^ (2.24) 



which proves the result flL13l) . 



It should be noted in conclusion that obtained solution is valid only for the considered expo- 



nential type distribution, eq. (|1.3|) . It is crucial for the above proof that the form of the distri- 
bution of a random element Jij doesn't change after its shift on a constant value. On the other 
hand, it is clear from the above proof that in the thermodynamic limit N — > oo the leading (in 
1/iV) contribution to is defined only by the very beginning of the distribution, P[J — > 0]. 
Therefore, the result En->oo = C(2) must be correct, also for the "rectangular" type distribution: 
P[0 < J < 1] = 1; P[J > 1] = (it is actually the model with this type of distribution which was 
studied in the replica solution U). For the discussion of other types of the matching models see 
e.g. 0] and references therein. 



3 Appendix A 



In this Appendix we prove that value of SE^ defined in eq. ( [2.19| ) is given by the rectangular N xk 
problem fl2~20l) . 

First, let us consider the most simple case k = 2: 



= E N - E N 



J 

(3-1) 

The above two problems, E^ and E^\ differ only by the permutation of two elements: (2,1) and 
(2,2), while all the other matrix elements in both matrices are the same. Nevertheless, even 
this "tiny" permutation, in general, can make the ground state configurations of the matrix S in 
the two problems to be quite different. Note that for the calculation of the above average energy 
difference 5E N we can make averaging over both simultaneously (keeping to be the same in 

~ (2*) (2) 

both problems) as well as separately for E N and for E N . 

For further proof it is important to introduce the concept of equivalence among the columns 
(and among the lines). We call the two columns j\ an d 32 ( or the two lines i\ and %<i) equivalent 
if the probabilities of the positions and (i, j 2 ) (or and (z 2 , j)) to be occupied in the 

ground state are equal. 

Due to obvious symmetry properties of the systems under consideration, it is evident that in 
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each of the above problems, E N and E N , all the columns on the right of the double vertical line, 
and all the lines below the double horizontal line are equivalent among themselves. On the other 
hand, the first two lines in each of the above problems are also equivalent between themselves, but 
they are not equivalent to the rest of the (N — 2) lines. Besides, in the problem E N we have first 
two columns which are equivalent between themselves, but which are not equivalent to the rest of 
the (N — 2) columns. Finally in the problem E N the first column is not equivalent to the rest of 
the (N — 1) columns. 

(2) 

The energy difference 5E N , eq(|3.1|), can be represented as follows: 



SE^= [^[JiM ~ Jm]j (3-2) 

where Jj^j and Ji(j)j represent the elements of the j-the column occupied in the ground states of 
the problems Eff and E$ correspondingly. 

In any ground state in each of these problems there are two columns in which the occupied 
elements belongs to one the first two lines. Therefore, in the summation over the columns in eq.( |3.2| ) 
we can find two, three or four columns in which one of the elements (or both) of the difference 
{Ji(jy — Ji(j)j) belongs to one of the first two lines. In the rest of the columns (their number can 
be (N — 2), (N — 3) or (N — 4)) both the element and the element belong to 

the rest of (N — 2) equivalent lines. Now, if with some probability we have the ground states in 
which the element is occupied in the problem E$ and the element is occupied in 

(2) 

the problem E N , then with the same probability we must have the inverse situation: the element 

~ (2) ~ (2) 

is occupied in the problem E N , while the element is occupied in the problem E N . 

Since both problems are defined by the same matrix (except for the two elements (2,1) and 
(2,2)) in the averaging over J^j in eq. ( |3.2| ) the contributions from all the columns with no elements 
in one of the first two lines must cancel. 

(2) 

Thus, to compute the difference d~E N , we should take care only of those columns which in the 

~ (2) (2) 

ground states of the problems E N and E N contain elements of the first two lines. 

Due to the equivalence of the first two lines and due to equivalence of the (N — 2) columns 
(j = 3, ...N) we can reduce all the relevant ground states of the problem E N to the following four 
non-equivalent basic configurations: 
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(a) 










(b) 




> 






(c) 










d) 







(3.3) 



where " •" represent the elements occupied in the ground state configuration of the matrix S, and 
"O" denote occupied element "0". Note that each of the above configurations represents the whole 
set of equivalent configurations. For instance, ( |3.3| (b)) represents all configurations with "•" in any 
of (TV — 2) positions (2, j), (j = 3, N), as well as all configurations with "•" in the position (2,1) 
and another "•" in any of (N — 2) positions (1, j), (j = 3, ...,N). The diagram ( |3.3| (c)) represents 
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all configurations in which one of the zeros is occupied. Note also that the configurations of the 
type: 






• 








• 






















(3.4) 



must be excluded from the consideration since they can not be the ground state being always higher 
in energy than the states represented in (fOl (d)). 

Due to equivalence of (N — 1) columns (j = 2, N) in the problem E$ here we have only two 
non-equivalent basic configurations: 
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* 












(a) 




> 






(b) 







Here for the occupied positions we use the notation 

A 2 ' 

J N 



(3.5) 



instead of "•" to distinguish them from 



(2) 

the ones in the ground states of the problem Ep 

Now to compute the contribution to the difference of the energies SE ( N ', eq. (|3.2|) , we have to 

(2) 

consider all possible combinations of the ground state configurations of the problem E N , eq. 

(2) 



and of the ones of the problem i% , ec L-( 

It is evident that if in the problem E]y we have one of the configurations of the type (|3.3| (a)) 

n(2) 



or (|3.3| (b)) and in the problem E N we have one of the configurations of the type ( |3.5| (a)) (all those 
in which no one "0" is occupied), then (since the two problems contain the same set of Jj/s) the 

(2) 



must coincide. Therefore, these two cases give no contribution to 5E 



N ) 



positions of " 
eq.®. 

It is also evident that the combination of one of the ground states of the type (|3T3"l(a)) or ( |3.3j (b)) 
with (|3.5|(b)) is impossible. For example, let us suppose that the ground state of the problem E N 
is the configuration ( |3.3| (a)), and the one of the problem E N is the configuration (|3~5|(b)). Then, 
according to the definition of the ground state, the energy of (|3.3| (a)) must be smaller than the one 

~ (2\ 

of the configuration ( |3.3| (d)), which in turn (since the problem E N contain the same set of Ji/s) 
must be smaller than the energy of the configuration ( |3.5| (b)). On the other hand, the energy of 



the configuration (3.3(a)) is equal to the one 










* 









* 

















(3.6) 



of the problem E N . Thus, the energy of ( |3.6| ) is smaller than the one of ( |3.5| (b)), and therefore 
(|3.5| (b)) can not be the ground state. 

Similar arguments shows that the combinations of ( |3.3j (c)) with ( |3.5| (a)), as well as ( |3.3| (d)) with 
(|3.5|(a)) are also impossible. 
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The combination of (|3.3| (c)) and ( |3.5| (b)) is allowed, but in this case, according to the definition of 
the ground state, the position of "•" in (|3.3| (c)) of the problem E$ must coincide with the position 
of "*" in Ob)) of the problem E {2) 
8Ef. 



N , and therefore this combination also gives no contribution to 



Finally we are left with the combination of the ground state configurations of the types (|3.3j (d)) 
and ( |3.5| (b)) which indeed give finite contribution to 5E^\ eq.( p.2|) . According to the discussion 



below eq. 



here we get no contribution from all the elements of all (N — 2) lines 



■AO, 



and thus all these elements fall out of the computation. The difference of energies between (|3.5| (b)) 
and ( |3.3| (d)) which does enter into the computation is equal to one of 2(N — 1) elements of the first 
two lines of the problem E N (note again, that the diagram (|3.5|(b)) represents one of 2(iV — 1) 

(2) 

configurations with one occupied "0"). Therefore, the average difference SE K N ' can be represented 
in the form of the "truncated" problem: 



SEP 
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(3.7) 



where the ground state is just one the configurations of the type ( |3.5| (b)). 

The generalization of the proof for arbitrary k is straightforward. First of all, it is evident that 
to get the proof for a general value of k it is actually sufficient to consider the problem with k = 3: 
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= E N - E N 



J 



(3-8) 

because here the first line and the first columns can "represent" all the other equivalent lines and 
columns of the problems with k > 3. 

Then we should make quite similar to eqs.( |3.3| ) and ( |3.5| ) classification (which turns out to be 
only slightly more cumbersome) of all non-equivalent ground state configurations of the problems 



E N and E^' according to the positions of the occupied elements of the first three lines. Simple 



,(3) 



analysis shows that here again the only relevant (for SE^ ) configurations of the problem E^' are 
the ones with all three zeros occupied, while in the problem E N these are the configurations with 
one or two of the zeros occupied. On the other hand, due to equivalence of the rest of (N — 3) lines 
(i = 4, ...N) one finds that all the elements of these lines fall out of the computation. In this 
way one obtains that energy difference 8E N is given by two or one of the elements of the three 

~ (3) 

first lines of the problem E N , which is just defined in terms of the "truncated" problem: 
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(3.9) 



The calculation of the actual values of 5E$ is presented in the next Appendix. 



4 Appendix B 



In this Appendix we prove that 



5E 



(fc) 
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k(k 



k-l 
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i=i 
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N-l 



(4.1) 



The solution of the case k = 2, eg. ( |3.7| ), is trivial. Here the ground state configuration is of the 
type ( p. 5| (b)), where the position of "*" must be at the smallest element out of 2(N — 1) non-zero 
Jy's. According to the distribution ( |1.3| ), for the average value of this element we get: 



6E 



(2) 
N 



1 



2(N 



(4.2) 



Now let us consider slightly more complicated case k = 3, eq. ( |3.9| ). Simple analysis of the the 
structures of possible ground state configurations shows that all of them can be taken into account 
in terms of only 3x3 matrix: 
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z 
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y 





X 



(4.3) 



Here x is the smallest element out of 2(iV — 2) equivalent elements J 2 j (j = 3,...,N) and J 3 j 
(j = 3, ...,N) of the second and the third lines; y is the smallest element out of two equivalent 
elements J21 and J31; z is the smallest element out of (N — 2) equivalent elements J\j (j = 3, ...,N) 
of the first line; the symbol "(g>" denotes the elements which do not enter into any ground state 
configuration. One can easily check that the matrix in eq. ( |3.9| ) can have only two ground state 
energies equal to x or equal to (y + z). Note that if we consider this problem in terms of the 3x3 
matrix ( f4.3|) , the element x could as well be placed in the position (3,2) (instead of "®" which then 
should be placed at the position (3,3)), as well as y could be interchanged with "(g)" in the positions 
(2,1) and (3,1). 



Now one can easily note that original 3x3 problem (4.3) is actually equivalent to the 2x2 
problem: 



E 






z 


y 


X 



(4.4) 
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where, according to the definitions of the random parameters x, y and z their statistical distributions 
are: 



P(x) = 2(N - 2) exp [-2(N - 2)x] 



(4.5) 



P(y) = 2exp(-2y) 



(4.6) 



p( z ) = (N - 2) exp [—(N - 2)z] 



(4.7) 



Keeping in mind further generalization of the solution for an arbitrary k, we solve the problem 
( |4.4j ) in the following way. Similarly to the procedure described in the beginning of Section 2, we 
can " shift" the elements of the second line (x and y) by the value of the smallest of them, and then 
integrate it out: 



+ 



1 



E 






z 





X 



(4- 



2(N~1) (N - 

The factor 1/(N — 1) in the second term of the above equation is the probability that y is smaller 
than x (if the smallest element is x, then the remaining problem will have all zeros at the diagonal, 
and the minimum energy of this problem is identically equal to zero). The solution of the remaining 
2x2 problem is trivial, and eventually we get the following result: 



6E 



(3) 
N 



3-2 



1 



+ 



N-l N-2 



(4.9) 



2(JV-1) 3(JV-l)(jV-2; 
Now the generalization of the above procedure for an arbitrary value of k becomes evident. First 



we note that all possible ground state configurations of the N x k problem SE^\ eq. (|2.20 ), can be 
taken into account in terms of the k x k matrix: 














Z\ 













Z2 
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Z(k-2) 


<g> 


® 




® 
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V\ 


V2 




V(k-2) 





X 



(4.10) 



Here x is the smallest element out of 2(N — k + 1) equivalent elements J(k-i)j (j = k,...,N) and Jkj 
(j = k, N) of the last two lines; yj (j = 1, (k — 2) is the smallest element out of two equivalent 
elements J(k-i)j and Jkj] Z{ (i = 1, (k — 2) is the smallest element out of (N — k + 1) equivalent 
elements Jij (j = k, ...,N) of the i-th line; and again the symbol "®" denotes the elements which 
do not enter into any ground state configuration. 

According to the above definitions of the random parameters x, {yj} and {z^} their probability 
distribution functions are: 
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P(x) = 2(N - k + 1) exp [-2(N — k+ l)x] 



(4.11) 



P(y j ) = 2exp(-2y j 



(4.12) 



P( Zi ) = (JV - k + 1) exp [-(JV - k + l) Zi \ (4.13) 
In this way we can reduce the calculation of 5E^ to the (k — 1) x (k — 1) matrix problem: 
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(4.14) 



/ 



Taking into account the equivalence of the first (k — 2) columns here we can integrate out the 
smallest element of the last line to get: 
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(4.15) 



/ 



Now one can easily see that all possible ground state configurations in the remaining (k-l)x(k-l) 
problem can be taken into account in the same way as in the previous k x k one, eq. (|4.10|) . Here 
we can reduce the number of relevant elements by choosing the smallest one between x and z^-2), 
as well as between each yj of the last line and J<k-2)j (j = 1> (k — 3)) of the previous line. In this 
way we get: 
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(4.16) 



where the random elements x and {yj} according to their definitions are now described by the 
following distribution functions: 



P(x) = 3(JV - k + 1) exp [-3(JV -k+ l)x] 
p (Vj) = 3 exp (-3%) 



(4.17) 
(4.18) 
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while the distribution functions of z^s remain unchanged, eq.( [4.13|) . In this way we can reduce the 
calculation of SE^ to the (fc — 2) x [k — 2) matrix problem: 
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n 2(N-1) ' (JV-1) 
Here again we can integrate out the smallest element of the last line to get: 



\ 
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1 k-2 
+ 



2(JV-1) (N-l) 



1 fc-3 

+ — — -E 



3(N-2) (N-2)' 
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(4.19) 
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Z/(fc-4) 





X 



J 

(4.20) 

Continuing these iterations till the last trivial 2x2 problem we eventually obtain the following 
result: 



5E 



(fc) 



1 fc - 2 

+ 



2(JV-1) (iV-1) 



1 



fc-3 



3(N-2) (N-2) 



1 jfc-4 
+ 



4(iV-3) (JV-3) 



1 



k(N-k + l) 



After simple algebra the above expression can be easily reduced to the following form: 



SEP 



fc(Jfe-l) 



fc- 1 



JV-1 JV-2 



JV- fc + 1 



(4.21) 



(4.22) 



which proves eq.( |4.1| ) 
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